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ABSTRACT 

It has been proposed that the Galaxy might contain a population of cold 
clouds in numbers sufficient to account for a substantial fraction of the total 
mass of the Galaxy. These clouds would have masses of ~ lCr 3 M and sizes 
~ 10 AU. We consider here the lensing effects of such clouds on the light 
from background stars. A semianalytical formalism for calculation of the 
magnification event rate produced by such gaseous lensing is developed, taking 
into account the spatial distribution of the dark matter in the Galaxy, the 
velocity distribution of the lensing clouds and source stars, and motion of the 
observer. Event rates are calculated for the case of gaseous lensing of stars in 
the Large Magellanic Cloud and results are directly compared with the results 
of the search for gravitational microlensing events undertaken by the MACHO 
collaboration. The MACHO experiment strongly constrains the properties of 
the proposed molecular clouds, but does not completely rule them out. Future 
monitoring programs will either detect or more strongly constrain this proposed 
population. 

Subject headings: galaxies: halos — galaxies: ISM — Galaxy: halo - 
gravitational lensing — ISM: clouds 
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1. Introduction 

The nature of the dark matter dominating the mass of the Galaxy remains elusive. 
Although nonbaryonic forms of dark matter are one possible explanation, it is possible that 
much or most of the mass may be baryonic. Recent results from gravitational microlensing 
experiments, such as MACHO and EROS (Alcock et al. 2000; Lasserre et al. 2000), have 
put severe constraints on the abundances of compact baryonic objects in the Galaxy, - such 
as brown dwarfs, old white dwarfs, neutron stars, or black hole remnants - thus virtually 
removing stellar objects as candidates for being the major constituent of the dark matter 
(Alcock et al. 2000; Freese, Fields, & Graff 2000; Lasserre et al. 2000). 

It has recently been proposed that the dark matter in the Galaxy could consist of 
small, cold clouds of H 2 (Pfenniger, Combes, & Martinet 1994; De Paolis et al. 1995a, 
1995b, 1996; Gerhard & Silk 1996). Walker & Wardle (1998) pointed out that such clouds 
could be responsible for the extreme scattering events (ESEs; Fiedler et al. 1994), which 
would occur when the line of sight to an extragalactic point radio source is crossed by 
the photoionized cloud envelope. In a later paper, Wardle & Walker (1999) considered 
the thermal stability of such clouds and showed that these clouds could be stable against 
heating by cosmic rays. These clouds could have masses of the order of Jupiter's mass 
(1O~ 3 M ) and radii R c i « 10 AU. If they exist, these clouds could also naturally explain 
the core radii of galaxies (Walker 1999) as well as the 7-ray emission from the Galactic 
halo (Kalberla, Shchekinov, & Dettmar 1999; Sciama 2000). Cloud-cloud collisions would 
steadily resupply the Galactic disk with gas to sustain steady star formation. 

The proposed clouds are not compact enough to produce any gravitational microlensing 
(Henriksen & Widrow 1995) but Draine (1998) demonstrated that they can produce 
magnification of background stellar sources. Refraction of light passing through the 
clouds would cause amplification of background stars in a way resembling gravitational 



-4- 



microlensing, which provides us with the possibility of using searches for gravitational 
microlensing events to either detect such gas clouds or constrain the properties of the cloud 
population. 

In this paper we calculate the lensing rate of the stars in the Large Magellanic Cloud 
(LMC) by clouds of molecular hydrogen, taking into account the spatial distribution of the 
clouds and all the relevant motions which contribute to this rate: motion of the observer on 
the Earth, proper motion of the LMC, the velocity distribution of stars in the Cloud, and 
the velocity distribution of lensing clouds in the Galactic halo. We suppose that the dark 
matter is composed of clouds of only one size and mass and that the clouds are transparent. 

The paper is organized as follows: in §|] we review some aspects of the lensing by 
gaseous clouds. In §|3| we derive basic formulae for the rate of lensing events produced by 
the gaseous clouds and for timescale distribution of this sort of lensing. In §|] we present our 
results in the form of parametric plots, covering a wide range of cloud models, and compare 
our results with those obtained by collaborations undertaking searches for gravitational 
microlensing events. We determine the region of parameter space which is not excluded by 
these experiments and other constraints. Finally, in §^ we compare our results with the 
limits placed on cloud models by other authors from different arguments. 

2. Physics of gaseous lensing 

The physics of lensing by gaseous clouds has been discussed by Draine (1998) and we 
repeat here only the salient points. The refractive index m is related to the gas density p by 

m(A) = l + a(A)p ; (1) 

a(4400A) = 1.243 cm 3 g- 1 and a(6700A) = 1.214 cm 3 g" 1 for H 2 /He gas with 24% He by 
mass (AIP Handbook 1972). 
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For small deflections, a light ray with impact parameter b will be deflected through an 
angle 

Let D s \ and D Q \ be the distance from source to lens, and from observer to lens (see 
Figure p]). If 60 is the distance of the lens center from the straight line from source to 
observer, then the apparent distance b of the image from the lens is given by the lensing 
equation 

b-b = D<p(b) , D = J^- . (3) 

For a point source, the image magnification is given by 

M(b) = M * , (4) 

6 1 - D(p'(b) 



CLU 







b 2 d 2 p z 2 dp 



(5) 



r 2 dr 2 r 3 cfr 

where r 2 = b 2 + z 2 . For a given bo there will be an odd number N(bo) of solutions hipo), 
i = l,...,N. The total amplification A(b ) = Y^i=iM(bj). The "trajectory" of the lens 
relative to the source is characterized by a "source impact parameter" p and a displacement 
x along the trajectory; for any x we have bo = (p 2 + x 2 ) 1 ^ 2 , and the "light curve" is just 
A{po) vs. x. 

Following Draine (1998) we define a dimensionless "strength" parameter 

S= ^ = 0.355 PMf^yf^ , (6) 
where M d ,R c[ , and (p) = 3M d /4:iiR d are the cloud mass, radius, and mean density 

As the pressure density-relation in such clouds is uncertain, Draine (1998) considered 
polytropic equations of state for polytropic index 1.5 < n < 5. (For n = 1.5 the cloud is 
isentropic, while for n — > 5 the central density becomes infinite). 
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For each polytropic index n there exists a specific value of S = S cr such that for 
S < S cr equation (|3|) has only one solution and magnification of the lens is always finite, 
even for b = 0. For S > S cr equation (fj) can have three solutions for sufficiently small 60 
and caustic lensing can occur. In this caustic regime the magnification becomes infinite 
for bo = and also for some finite b$ c given by the condition that at this b$ c two of the 
roots of equation coincide. At b = b 0c light curves exhibit conspicuous caustics with 
magnification going to infinity. Cases of caustic and noncaustic lensing are illustrated in 
Figure |2] for n = 1.5. 

For S < S cr , we can obtain a simple analytical formula describing the dependence of 
the central magnification M(0) upon the strength parameter S. If we define dimensionless 
variables p = pj (p) and f = r/R d , it is straightforward to show that as b — » 

b - k ^-2S b j d 4§ (7) 
J r ar 


The integral in this expression depends only on the cloud density profile. 
Substituting this result into (f|) we get 

M(0) = (8) 

(1 -s/ S cr f 



where 



1 



-1 







Thus M(0) — ► 00 as S — > S cr , but it is finite for smaller S. At S = S cr the number 
of solutions of equation (|3|) changes from 1 to 3 (for b = 0) and for S > S cr we enter the 
caustic regime. For n = 1.5 S cr = 0.026 and for n = 4.5 S cr = 1.72 x 10~ 6 . 



-7- 



3. Rate of gaseous lensing events 

When a cloud crosses the line of sight between the distant star and the observer on the 
Earth, it amplifies or attenuates the brightness of the star [in fact each gaseous lensing curve 
exhibits some demagnification before and after the peak magnification (Draine 1998), but 
we will be primarily interested in amplification in this paper]. This effect can be observed 
by searches for gravitational microlensing in our Galaxy which monitor the brightness of a 
large number of stars. 

Let D s be the distance to the source stars being lensed by gaseous clouds in the Galaxy. 
Consider a lens located at some distance Dol = xD$ from the observer; the distance from 
the lens to the source star is D LS = (1 — x)D s , and D = x(l — x)D s . 

Both lensing regimes, caustic and noncaustic, are important. Since the lensing 
parameter 

maximum S is attained when the lens is placed midway between the source and the 
observer, and is equal to S max = oc(p)Ds/(4:R c i). As x changes towards x — or x — 1, S 
declines to zero. Thus, if S max > S cr , there are critical values of x 



such that for X\ < x < x 2 lensing is caustic while for x < x x and x > x 2 lensing is in the 
noncaustic regime. If S max < S cr lensing is always noncaustic. 

For the observer to see a noncaustic lensing event with a magnification larger than 
some threshold magnification A t , the lens has to pass near the source star with a sufficiently 
small (unlensed) impact parameter 6 - in other words, magnification M > A t if b < b ot (A t ), 
where bot(A t ) is given by 

M(bot(A t )) = A t . (12) 
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Fig. 1. — Geometry of lensing. S is the source, the center of the lens is at L, and O is the 
observer. 
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As the cloud moves through the sky, all the source stars lying in a strip on the sky 
along the lens trajectory with angular width 25 given by 

5(A t ,x)=b ot (A t )/D OL , (13) 

are amplified with magnification M > A t . 

Let vs, Vx,, and v G be the transverse velocities (i.e. velocities perpendicular to the 
line from observer to source star) of the source star, the lens, and the observer. Then 
the relative transverse velocity vj_ of the gaseous lens and the source star as seen by the 
observer is just 

v± = v L - [xv s + (1 - x)v ] (14) 

(Han & Gould 1996). 

It is important to distinguish events with different durations. Any real gravitational 
microlensing experiment has some finite probability < 1 of detecting a lensing event which 
occurred during the monitoring campaign and this detection efficiency function depends 
strongly upon the timescale of the observed event (Alcock 1997). We define the timescale r 
of the lensing event as the time which the light curve spends above a threshold magnification 
A t : 

T = M ' (15) 

where ip is the angular impact parameter of the lens's trajectory relative to the source star; 
(p < 5 (otherwise M < A t ). 

In calculations of the rate of lensing we must take into account the fact that lensing 
clouds and source stars have some distribution in velocity space; we average the lensing rate 
over these distributions to obtain the true expected rate. The general formula for the event 
rate per lensing cloud is 

8 

dN = 2S 5 J d<p J rfv L J dv s / L (v L )/ s (v s )^0 T (r(y9, |vj)), (16) 
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Fig. 2. — Plots, showing lightcurves for n = 1.5 with different ratios S/S cr : (from top to 
bottom) 4, 2, 1.2, 0.4, 0.12. On the three top panels one can easily see caustics produced 
by multiple images. For lightcurves corresponding to other values of n and S/ S cr see Draine 
(1998). 
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where fa and f$ are the velocity distribution functions for lensing clouds and source stars, 
respectively, is the number of source stars per unit solid angle on the sky, and T is the 
detection efficiency function for a lensing event with timescale r. 

We assume for simplicity isotropic Maxwellian distribution functions for the velocities 
of both lenses and source stars with dispersions ol and as respectively, taking into account 
the fact that the object containing stars is moving in general, i.e. there is some offset 
velocity v c in the latter distribution: 



1 



2 l„2 



h(v L ) = —^e-^, (17) 
and 

fsM = ^7^ e " (VS " Vc)2/4 - ( 18 ) 
Our calculations are not very sensitive to this assumption in the sense that the order of 
magnitude result does not depend strongly upon the exact shape of the velocity distribution. 
The important assumption of isotropy of the distribution function permits analytical 
simplifications. 

Consider a coordinate system with the z-axis lying along the line of sight from observer 
to the source, x-axis perpendicular to the line of sight so that observer's velocity vq lies in 
the x2-plane, and y-axis perpendicular to those two, so that vo y = 0. Taking distributions 
(|T7|) and flIB|) it is shown in Appendix |A1 that the total event rate per source is 

Ntot = 4£>* j dx XnLi f[ ix [ A ^ ch e- c Ju\(uv ch )e~ Au2 I (Bu) du, (19) 

L b 

where Jo is the modified Bessel function of order zero, nz(x) is the number density of 
lensing clouds, v c h(x) is some characteristic velocity (the final result does not depend upon 
the choice of this velocity), function k is defined as 



(20) 
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and A(x), B(x) and C(x) are given by 

A = 2 V < CH 2 2 » # = 2^2 2 + (^CT + (1 ~ S^Os) 2 , 

erf + aHcr^ erf + x^crg v y 

n _ x 2 v% + (xv cx + (l-x)v 0x ) 2 

°~ ^TTx^l • (21) 

This formula for the lensing rate is quite general. For example, it is directly applicable 
to the case of gravitational microlensing [with corresponding 5(A t ,x)}, in which case it is 
identical to the formula for the event rate obtained by Griest (1991). 

It is interesting to know the timescale distribution of the event rate to have some idea 
of where to look for the events produced by the gaseous lensing. The timescale distribution 
is given by (see Appendix |A|) 

^ = 32D s ^i fd ^f^'f'K -c f^L=e-^I„(B«)du, (22) 
dr t J a i + x 2 ai J Jl-u 2 

o ^ 15 o v 

where b ot (x, A t ) is defined in fll2|) and v c h, entering the definitions of A, B, and C in (pT[) , is 
given by v c h = 2b 0t (x,A t )/r [compare with a similar formula by Griest (1991) for the case 
of gravitational microlensing]. 

Closely related to the timescale distribution dN tot /dr is its integral from r min to 
infinity, that is the frequency of events with timescales exceeding r„ 



1 mini 

oo 



A",„,(r > -„„„) = / ^-^dr 

dr 



32^ / ^ M*)^A) e -c } u 2~ {u)e -Au> h {Bu) ^ (23) 



' mm q 



where function R is now 



k( u )= I (j) T (^VT^)ds. (24) 



If one wants to know the true event rate, for 100% efficiency, one need only set <f> T = 1 



in formulae (|22"D and (|23|), as we will do in our consideration of the event rate distributions. 
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In the next section we show typical plots of N tot {r > r min ) for some clouds characterized 
by polytropes with n = 1.5 and n = 4.5. 

4. Results 

In this section we consider the lensing of stars located in the LMC and compare with 
the results obtained by the MACHO collaboration, which has monitored 11.9 million LMC 
stars for 5.7 yrs (Alcock et al. 2000). 

4.1. Parameters of the event rate calculations 

The LMC has a distance D LMC ~ 55 kpc from the Sun, and its position in galactic 
coordinates is b = —32.8° and I = 281°. Jones et al. (1994) finds the motion of the LMC to 
be 

vi = —16 ± 60 km s _1 

v b = 328 ± 60 km s" 1 (25) 
v rad = 250 ± 10 km s _1 . 

We use these proper motion measurements for obtaining our event rates. The velocity 
dispersion of the lensing clouds was assumed to be <tl = 220 km s _1 and that of the stars 
in the LMC was taken to be a s = 10 km s -1 . 

We performed our calculations both for the case of the true event rate (for 100% 
detection efficiency function) and also for the the case of the MACHO detection efficiency 
function which is given in Alcock et al (1997). In the first case we assume a threshold 
magnification A t = 1.1 and in the second we take A t = 1.35, which is the current threshold 
magnification of the MACHO experiment (Alcock et al. 2000). Estimation of the MACHO 
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detection efficiency function T is problematic because it is known in terms of the duration 



of a gravitational microlensing event, which differs significantly from our definition (15) 
for the gaseous lensing timescale. The method by which we estimate <p T is described in 
Appendix ||, and our adopted (fi T is shown on Figures [| and f|. 

We model the lensing clouds by self-gravitating H 2 -He (nonrotating) polytropes of 
radius R c [. The polytropic index n (T oc p 1 /™) must be in the range 1.5 < n < 5; for n < 1.5 
the cloud would be convectively unstable, while for n > 5 the central density is infinite. 
We do not expect T(r) and p(r) to be accurately described by a polytropic model, but a 
slight rise in temperature toward the interior may be reasonable since the interior will be 
heated by high energy cosmic rays, with the few cooling lines [e.g., H2 0-0S(0) 28.28pm, 
or HD O-OR(l) 112 pm] very optically thick. Furthermore, polytropes have T(R C {) = 0, so 
the density structure near the surface is unphysical. Table 1 of Draine (1998) gives various 
properties for H2-He polytropes. In our current calculations we consider gaseous clouds 
with polytropic indices n = 1.5 and n = 4.5, to compare extremes of behavior^. 

Following Widrow & Dubinski (1998) we adopt a Navarro, Frenk, & White (1996) 
model of a spherical halo, composed of gaseous clouds: 

" i(r) = ! U«M d r{r + a,y (26) 

where r is the distance between the galactic center and the point of interest, a s = 15.9 
kpc is the core radius, and < < 1 is the fraction of the total mass of the Milky Way 
M MW = 6.5 x 10 11 M contributed by the population of self-gravitating H 2 clouds. 



2 We assume that it is the density-radius relation of the clouds which is determined by the 
polytropic law, not the behavior of the gas on dynamical timescales (otherwise the n = 4.5 
polytrope would be dynamically unstable.) 
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Table 1. Models used for timescale distribution. 



model 


n 


Md 1 


Rci 1 


rp 2 
1 c 


2GM cl m H /3k B R cl 3 


Smax 




A 


1.5 


10~ 2 


6.6 


200 


109 


2.6 x 10" 


Q 

-o 


B 


1.5 


10" 2 


2.6 


507 


276 


0.106 




C 


1.5 


10" 2 


1.1 


1200 


652 


3.3 




D 


1.5 


10 -3 


3.7 


36 


19 


2.6 x 10" 


-3 


E 


1.5 


10 -3 


1.5 


88 


48 


0.096 




F 


1.5 


10 -3 


0.6 


220 


119 


3.8 






l.O 


1U 


n s 

U.O 


1 7 


Q 


n 1 9 




H 


1.5 


io- 4 


0.4 


33 


18 


1.9 




I 


4.5 


lO" 2 


64 


13 


11 


2.9 x 10" 


-7 


J 


4.5 


io- 2 


8.8 


93 


81 


8.1 x 10" 


-4 


K 


4.5 


io- 2 


1.3 


628 


552 


1.7 




L 


4.5 


io- 3 


5 


16 


14 


7.8 x 10- 


-4 


M 


4.5 


IO" 3 


0.7 


117 


102 


2.0 




N 


4.5 


1(T 4 


0.4 


20 


18 


1.9 





1 M d is in M Q , R d is in AU. 

2 Cloud central temperature, in Kelvins 

3 Temperature at which the thermal energy of a hydrogen atom at 
the cloud surface equals its gravitational binding energy, in Kelvins 
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4.2. Timescale distributions 

On Figures |3| and [| we plot the dependence of the integrated timescale distribution 
Ntot( T > T min), as a function of T min , for some representative values of the cloud mass M d 
and radius R c i. The cloud parameters are given in Table [3] and we assume <ft T = 1. 

All timescale distribution curves exhibit a plateau for small timescales and then a rapid 
falloff at larger event durations. From these plots we see that the characteristic timescale 
of the lensing events is quite small, varying from several days to several tens of days. This 
is as expected, because cloud sizes are small, ~ 10 AU, sufficient amplification occurs only 
when 6 < (0.1 — 0.3)R c i, and for typical transverse velocity ~ 200 km s _1 we get just this 
range of event durations. 

The time distribution seen by a real experiment would be modified by the detection 
efficiency, which may be small for short timescale events. In addition, events lasting > 50 
days, which occur when the transverse velocity of the lens is quite small, will be affected 
by changes in the velocity of the Earth over this period (the so-called "parallax effect" 
(Paczyhski 1996; Gould & Andronov 1999)). 

4.3. Results for the event rate 

Calculation of the event rate is quite straightforward in the case of noncaustic lensing 
because in this case foreground star crossing the cloud produces magnifications A > A t 
in one continuous span (see Figure |2| for the cases S < S cr ). In the caustic regime the 
situation is more complicated, since for clouds located at the middistance between observer 
and foreground star these magnifications A > A t will be not only in the region of the closest 
approach but also at the caustic spikes (see Figure || for the cases S > S cr ) and we want to 
account for this. 
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0.1 1 10 



r , days 



Fig. 3. — Distribution of r for lensing by selected n = 1.5 polytrope clouds. The timescale 
r is defined to be the time during which the magnification A > A% — 1.1. The rate of events 
with timescales larger than tq is shown, for detection efficiency function <p T = 1. Timescales 
of a few days are typical; for model F, for example, 50% of events have r < 3 days. The 
solid line labelled <j) T shows the detection efficiency function used in the MACHO event rate 
calculations in Figures [7| and ||. 
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Fig. 4 - 



Same as Figure |3| but for n = 4.5. 
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Fig. 5. — Contour plot of the rate of events with peak magnification A > At — 1.1 and any r 
produced by gaseous lensing as function of M d and S max = 3aM d D LMC /16nR d . Each cloud 
is assumed to be an n = 1.5 polytrope. Thick solid lines represent levels of constant event 
rate labelled in units of 10~ 6 events yr" 1 per source, for f d = 1. Dash-dotted lines are the 
lines of constant R d in these coordinates, labelled by R d in AU. To the left from the thick 
dashed line "rate = 0" the clouds are unable to produce magnification A max > A t = 1.1 
(see equation fl2"8|)). Dashed vertical line at S max = S cr separates caustic and noncaustic 
regimes. Dotted lines show the lower limit M c \ / R c i so that a molecule at the surface will be 
bound for an atmospheric temperature of 10 K or 15 K. Dashed line at M = 0.02 M is the 
upper limit on M d suggested by Wardle & Walker (1999). The unshaded region is allowed. 
Squares with letters show the positions of the sample cloud models from Table |l[ Note the 
very high event rate for almost all cloud models (to the right from line "rate = 0"). 
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We adopted a simple approach for the caustic regime. One can see from Figure |2| 
(upper panel) that for a given polytropic index there exists a specific value of S* > S cr such 
that the inner minimum of the lightcurve (between the central part and one of the caustic 
spikes) has magnification A t . For example, in the case n = 1.5 S± = 0.112 for A t = 1.1 
or S* = 0.101 for A t = 1.35. It is obvious that when S cr < S < S+ the picture is similar 
to the noncaustic case since magnifications A > A t are reached in one continuous span. 
If S > S* there are also caustic spikes producing events but one can easily see that the 
duration of events caused by these spikes decreases rapidly with growing S (Draine 1998) 
and we simply ignore them in our calculations. That is we assume for S > S± that A > A t 
only in the central part of the lightcurve. This gives us only a lower limit on the event rate 
but a useful one since the number of events produced by caustics is relatively small and the 
caustic spikes are of short duration. 

In Figures [5] and |6| we show the true event rate (i.e. for efficiency T = 1) due 
to gaseous lensing for various cloud parameters and for two polytropic indices of the 
clouds: n = 1.5 and and n = 4.5. We characterize the clouds by mass M c \ and S max , 
the "strength" parameter for a cloud located midway between the Earth and the LMC: 
Smax = oi{p)Dlmc /^Rci, where Dlmc = 55 kpc. For each value of S max and M c i the cloud 
radius R c i is obtained from equation (^j): 

«^f 3aM f tMC ) V4 = 1.25x10" cm 1 (27) 

Lines of constant cloud radius are plotted on each of Figures ^|-^| as dot-dashed lines. We 
consider here events with all timescales and magnifications larger than A t = 1.1. 

Our adopted bound on M r j comes from two considerations. First, clouds which are too 
massive would be similar to the objects known to be unstable to collapse to form low mass 
stars. Second, Wardle & Walker (1999) obtained an upper bound on M r j < 1O~ 1 ' 7 M by 
considering a cooling mechanism which could stably radiate away the heat deposited by 
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cosmic rays. 

We restricted the parameter space of S max and M c i because of the following reasons. 
The upper bound in S max is just 10 because this corresponds to a cloud with M c i = 1CT L7 M 
and radius R c i = 1 AU. We do not consider smaller radii because extreme scattering events 
favor clouds with sizes of the order of several AU (Draine 1998). 

The minimum value of S for a cloud to be detectable by light curve monitoring with a 
threshold magnification A t can be obtained from equation @: 



where S cr is given by equation (§). For S max > Sdet the cloud population will produce a 
nonzero rate of events with A max > A t . 

Dotted lines show the restriction imposed by the requirement that the gravitational 
binding energy of gas at the surface exceed the thermal energy for an estimated atmospheric 
temperature T (Draine 1998): 



Two boundaries are shown, for atmospheric temperatures T = 10 K and T = 15 K. The 
shaded region below these curves is prohibited. 

The thick solid lines are contours of constant event rate. The lensing rate is very high 
for all the cloud models appreciably to the right of the "rate = 0" line. 

Figures |7| and || show the results for the lensing rate for polytropes with n = 1.5 and 
n = 4.5 for a detection efficiency function approximately that of the MACHO experiment. 
All the notation is the same as in the case of <p T = 1. The "rate = 0" boundary has moved 
(relative to Figures [| and []) because we now require a threshold magnification A t = 1.35. 
Small glitches on the contours are artifacts of the numerical procedure. We consider here 




(28) 




(29) 
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Fig. 6. — Same as Figure |5| but for n = 4.5. 
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the events with all timescales filtered through the detection efficiency function discussed in 
Appendix ||. 

The rate which would be seen by the MACHO experiment is significantly smaller 
than the rate for A t = 1.1 and <p T — 1, sometimes by two orders of magnitude. This is 
primarily a consequence of the reduced sensitivity to short timescale events of the MACHO 
detection efficiency function. Nevertheless the predicted rate is still quite high and in the 
next subsection we compare the predicted rate to that actually observed. 

4.4. Comparison with the MACHO results 

How shall we compare the results of the MACHO experiment with our simulated plot of 
what it would see if the dark matter is composed of molecular clouds? First of all, we need 
to estimate what fraction of the events actually observed by MACHO could be attributed 
to gaseous lensing. This can be done by noticing that the typical timescale of gaseous 
lensing events is very small, while the shortest timescale event reported by MACHO has 
t = 34 days. For the types of timescale distribution seen in Figures ||] and £| it is extremely 
unlikely that there are more than 1 or 2 gaseous lensing events among those observed by 
MACHO, since otherwise a large number of short timescale events (with durations less 
than 10 days) would be seen, in conflict with the distribution of timescales observed by the 
MACHO and EROS (Alcock et al. 2000; Lasserre et al. 2000) collaborations. It would 
also be completely inconsistent with the EROS and MACHO combined limits on the rate 
of very short timescale events, from 15 minutes to several days, for which they claim that 
the analysis of two years of data on 8.6 million stars found no short-duration "spike" events 
(Alcock et al. 1998), implying an upper limit N < 10~ 7 events yr _1 . 

The MACHO collaboration observed 11.9 million stars in the LMC for 5.7 years. This 
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Fig. 7. — Contour plot of the event rate of gaseous lensing for various parameters of the 
lensing clouds, for the MACHO detection efficiency function. Polytrope n — 1.5 is assumed 
and threshold magnification A t = 1.35. All the symbols and lines are the same as in Figure 
|5j. Note the substantial decrease of the event rate due MACHO's reduced sensitivity to short 
timescale events. 
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means that the rate which could possibly be due to gaseous lensing events can be at most 
(0.01 — 0.03) x 10~ 6 events yr" 1 . It is clear from our simulated observations of gaseous 
lensing on Figures [7] and |8] that if j c \ > 0.1 the contour describing such a rate will lie very 
close to the line S max = Sdeti where the event rate drops to 0. Any cloud model with 
Smax < Sdet will of course be undetectable and thus be allowed. Models with S max > l.lSdet 
are prohibited, because otherwise MACHO would have seen an enormous number of gaseous 
lensing events. 



These considerations and the restrictions described in § |4.3| are combined together to 
produce an exclusion plot (Figure |9|). One can see that an allowed region exists where 
clouds basically cannot be seen by the MACHO experiment and are not prohibited by 
the maximum mass and "thermal evaporation" constraints. For larger n (softer equation 
of state) the allowed region shrinks. It is likely that for n — > 5 this region disappears 
completely, though we did not run our calculations for n larger than 4.5. 

Also in the case n = 4.5 another allowed region appears for small cloud radii (see 
Figure |8[). It is bounded by the level contour corresponding to the rate 3 x 10~ 8 events yr -1 
and goes into the region of the small cloud radii (see Figure | b). 

The allowed regions can be further narrowed by noting that clouds are highly unlikely 
to have high central temperatures. Detailed studies of their thermal structure are not 
available, but it seems reasonable to exclude models of clouds with central temperatures 
> (50 — 100) K. For n = 1.5 this places another limit within the allowed region (see Figure 
|a). In the case n = 4.5 the allowed region for small radii virtually disappears because of 
this restriction (see Figure ||b). 

In the absence of realistic models for the thermal structure of these clouds, we can 
only use n = 1.5 polytropes as the most conservative model from the standpoint of gaseous 
lensing. For this case, we see that the allowed region is quite substantial, and contains the 
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Fig. 8. — Same as Figure |7| but for n = 4.5. Notice another allowed region at small cloud 
radii, bounded by contour corresponding to the event rate 3 x 10~ 8 events yr -1 and line 
M = 1(T L7 M m . 
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model with M d = 1O~ 3 M and R d = 10 AU which was favored by Draine (1998). 

It is also possible to expand the allowed region by assuming that the clouds contribute 
only a small part f d <C 1 of the dark matter (see § f4.1|) , but to get any significant effect we 
need to suppose that this fraction is < 10~ 2 which would rule out these clouds as a main 
constituent of the dark matter, and eliminate them as the explanation for the extreme 
scattering events. 

It is clear that our lensing restriction represented on Figure || is sensitive to the adopted 
density profile of an individual cloud, but it is relatively insensitive to changes in the spatial 
and velocity distributions of the dark matter, since changes in these parameters are unlikely 
to produce variations in the event rate of > 2 orders of magnitude, which would be needed 
to move the restriction given by the observed lensing rate significantly away from the line 
"rate = 0". 

5. Discussion 

Wardle & Walker (1999) suggested that heating caused by cosmic rays in the cold 
molecular clouds can be balanced if particles of solid H 2 can form. They argue that these 
particles cool the cloud by thermal continuum radiation if the cloud mass lies in the 
range 10~ 7 ' 5 — 1O _L7 M . Our calculations show that the lensing event rate would grow to 
extremely high values for masses below 1O~ 4 M (see Figures [7| and |8| ). 

Kalberla, Shchekinov, & Dettmar (1999) proposed that the 7-ray background emission 
from the Galactic halo seen in EGRET (E > 100 MeV photons) data can be produced 
by interaction of high-energy cosmic rays with small dense clouds. They modeled this 
effect with different types of spatial distribution of dark matter in the form of these clouds 
and found that their best fit to the EGRET data for uniform density clouds occurred for 



-28- 




Fig. 9.— Exclusion plot for cloud radius R c i and mass M c i, showing region of allowed 
parameters (unshaded) and prohibited (shaded regions) for clouds with polytropic index 
a) n — 1.5 and b) n — 4.5. Dotted lines represent limitation set by requirement that 
the atmosphere be gravitationally bound. Dashed line shows mass restriction by Wardle & 
Walker (1999). Region above the dash-dotted lines contains models with central temperature 
larger than 50 K and 100 K, respectively, which might exclude them (line shading). Note that 
another allowed region appears for small cloud radii in the n = 4.5 case, bounded by dashed 
line roughly corresponding to the event rate 3 x 10~ 8 events yr _1 and line M = 10~ L7 M & . 
As discussed in the text, the n — 1.5 polytropic model is the most conservative assumption 
in the absence of detailed thermal models for these clouds. Star on the n — 1.5 plot shows 
the position of the cloud model with R d = 10 AU and M d = 1O~ 3 M . 
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Md/R^ ~ 10 3 M /(6 AU) 2 . Their conclusions about the cloud parameters are model 
dependent, but the quoted values fall within the allowed range on Figure |9] for n = 1.5. 

Though it is impossible to observe the gravitational microlensing by the molecular 
clouds in our Galaxy it becomes possible if one observes them in other galaxies. A search 
for such lenses in the Virgo cluster is now underway (Tadros, Warren, & Hewett 2000) and 
preliminary results indicate that galaxy halos are unlikely to primarily consist of object 
with masses smaller than 10~ 5 M Q . We confirm this restriction but also place more rigorous 
constraints on the cloud properties since we can reject many models with higher cloud 
masses. 

Another possible approach would be to look for periods of demagnification in the 
light curves of lensing events. For instance, the OGLE-II survey (Udalski et al. 1997) has 
collected a large sample of light curves for stars in the Galactic bulge, including ~ 300 
lensing events (Wozniak 2000). Of these, two or three appear to show demagnification 
preceding or following the central peak in brightness, although the statistical significance 
of the demagnification has not yet been established (Wozniak, private communication). 
Demagnification cannot be produced by gravitational lensing by one or more point masses, 
but would be a characteristic signature of gaseous lensing. 

6. Summary 

In this paper we have considered in detail the lensing of stars in the LMC by small 
self-gravitating molecular clouds which have been proposed as a candidate for the dark 
matter in the Galaxy. Lensing would occur because of the refraction of optical light in the 
clouds, with resulting magnification of the source. We have developed a semianalytical 
formalism for calculation of the lensing rate including the spatial distribution of the dark 
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matter, finite velocity dispersions of the lensing clouds and the source stars in LMC, and 
proper motion of the observer and LMC itself. 

Our calculations were carried out for a single mass and size cloud population. One can 
easily extend the analysis to a distribution of cloud masses and sizes by simply convolving 
event rates obtained for a uniform cloud population with the appropriate distributions. 

Lensing events might be detectable by searches for gravitational microlensing. This 
has allowed us to obtain constraints on the cloud properties by comparing our calculations 
with observational results obtained by the MACHO collaboration. 

We found that almost the only possibility for the dark matter to be in the form of such 
molecular clouds is for the clouds to be sufficiently weak lenses so that their lensing effects 
are below the detection threshold of the MACHO experiment, since otherwise a very large 
gaseous lensing event rate would have been detected. This still leaves an allowed region in 
parameter space where these clouds could exist and not contradict the limitations posed by 
lensing experiments and simple physical considerations. 

Though we have performed event rate calculations for only one halo model [given by 
equation (f^j)], we expect our results to be relatively insensitive to the particular form of 
the spatial distribution of these clouds, or the assumed shape of the velocity distribution of 
the clouds and the stars in the LMC. 

Future microlensing experiments with a lower detection threshold magnification A t 
could either detect these clouds or strengthen the constraints on their properties. 
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A. Derivation of equations ([19]) - (|23| ) 

For the purposes of the lensing rate calculations the velocities of the lens, observer, 
and source along the line of sight are not important, because only transverse motions are 
significant. It means that we can immediately perform integrations in equation (|16D over 
v Sz and v Lz . 

Taking into account ([15]) we can perform integration over ip in fllTf ) and define the 
result function k(v±): 

kM S Mt(<p))*P = }<Pr ^ X ^ DoL VT^"j ds, (Al) 

so that now 

oo oo oo oo 

dN = 25—^- / dv Lx / dv Ly / dv Sx / dv5»|vj.|/ L (vi, x , v Ly )f s (v Sx , v Sv )k(\v±\). (A2) 



'OL 

— oo — oo — oo 



D, 

We now change variables from vl x and vl v to v x and v y via equation ([i4|) : 

oo oo oo oo 

dN = 28 2 ff 2 a 2 D / dv * J dv V f dv Sx J dVsyy/vl + + V$) 

1 ^ L^ S OL 



-oo — oo — oo — oo 

,2 i /„, „,„, \2 /„, „, \2 i /„, „, \2' 



x ^ (Vx ~ (1 ~ ~ X^Sz) + K ~ XV Sy) _ ( v Sx ~ Vex) + ( v Sy ~ VcyY j_ 



We integrate over t>s x and i>s y to obtain after lengthy but straightforward calculations 
7r(a£ + x 2 a 2 s )D OL 



00 00 

rfA> = 2£ 5 / dva, / dv y Jv| + v2«(J^+«2) 



-00 —00 
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x ex P {- a 2 L X + x V 2 V a 2 + Ci«» + Ca«y " c} , (A4) 



where 

x 2 t> 2 „ + (xv cx + (l-x)v 0x 

o\ + I 2 (T S 



C = cy v " ' 2 y ^ , (A5) 

c 

cr^ + z: 2 0-| J " erf, + a^Og 



^ _ n XV cx + (1 - x)^ _ ^cy 



Introducing polar coordinates in velocity space v x = v cos (3, v y = v sin /3 and 
normalizing t> to some characteristic velocity v c h we get 

v 3 00 2?r 

c/A> = 25—, SV , ch ,^ e~ c [ du u 2 K(uv ch )e- Au2 [ e ™c h (b <** P) > d(3 ( A7 ) 

ir(ai + x 2 a 2 s )D L J ' V 



where 

a£ + x^ 

The last integral in equation (|A~?| ) can be reduced to 



A = TOT? ■ ( A8 ) 



2tt 2tt 

e 

o o 



WAy/g+GcosV+fio)^ = f e uv ch ^(lH!cosf3 d p = 27r/o ( Bu j ^ ( A g) 



where Jo is the modified Bessel function of order zero, and 

&i + x 2 a s 



B = + C2 2 = 2 ,^2 2 Jx 2 < + (x^ + (1 - x)^) 2 = 2VIC. (A10) 



So, finally we get for the event rate produced by one lens 

y 3 00 

c/A>(x) = 4,5— e~ c / «Vm; c/l )e- Au2 /o (Bu) du. (All) 

{ai + x 2 a I s )D L J 

Thus, the total event rate per source, produced by all lenses between us and LMC, is 
$iven by formula (|19|). 



Now, we obtain the formula for the timescale distribution (|2%p . Note that r and y> 
- the angular distance from the lens's trajectory to the source star in the perpendicular 
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direction - are directly related by equation (|15|): 



d,T u 2 

dip = -6(x, A t )— j. _ , (A12) 



T y/T 



u 



with u = v±/v cn , v c h = 2b Qt (x, A t )/r. To obtain dN tot /dr we simply omit the integration 
over ip in flT5|). Also, the velocity v± is limited by 2bot(x,A t )/r. All the other steps are 
analogous to those which we have done in derivation of (|19|) and we obtain equation (|22|). 



If we are interested only in events whose durations are larger than some chosen value 
T min , then to get the total rate we should integrate in ( Jl6| ) over ip not up to 8(x, A t ), but up 
to 



^ma, = \ 5 2 - = 5VT^, (A13) 

with u = v±t/25xD. One can easily see that we will get in this case for N tot (r > r min ) the 
formula ([23]). 



B. Detection efficiency function 

The detection efficiency function </> characterizes the sensitivity of the experiment to 
events with different durations. We will use the detection efficiency (j)f(A, i) assumed by 
MACHO experiment to estimate how it would affect gaseous lensing (see Figure 8 in Alcock 
et al. (1997)). 

We first recall some simple facts about gravitational microlensing. The amplification 
M in the case of gravitational microlensing is determined by the formula 

M(u) = " 2 + . 2 W2 , u = (u 2 mm + te 2 /?) 1 ' 2 



r E v± 



AGmDx{l -x) ll/2 



c 2 



(Bl) 



where b m i n is the impact parameter, m is the mass of the lensing object, t is the time of 
observation, and v± is, again, the relative velocity of the lens and the source star. 
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Recall that we defined the timescale r to be the time spent with magnification M > A t . 
From equation (|B1| ) we find that, for gravitational microlensing, 

1/2 



At A max 



= t 



2 _ 2 

U t U min 



(B2) 



J 



where Ut is the value of u at which the magnification is equal to A t and A max is the 
maximum magnification for a given u m i n . 

We use a very simple algorithm for converting detection efficiency function from t to 
r: since r oat, their ratio r/t is a function of A maa . (or w m j n ) for a given A t . Thus we can 
average this ratio over the distribution of the maximum magnifications A max and assume 
that the value of known 4>i(t) is equal to the r (< t >), that is 

Mr)=^i{r(<r/i>y i y (B3) 

The calculation of < r/t > is quite straightforward and can be most easily done in 
terms of u min (see the last equality in equation (|B2p). Since the distribution of lensing 
events over u min is just constant 

Ut 1/2 

/ K - <J du min ( * 



For A t = 1.35, accepted by MACHO experiment, < r/i >= 0.776, which means that 
r (r)=0 £ (1.288r). 
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